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We have carried out first-principles calculations of Be(0001) thin films to study the oscillatory 
quantum size effects exhibited in the surface energy, work function, and binding energy of the atomic 
hydrogen monolayer adsorption. The prominent enhancement of the surface density of states at the 
Fermi level makes Be(0001) thin films more metallic compared to the crystalline Be. As a result, the 
calculated energetics of Be films and the properties of atomic H adsorption onto Be(0001) surface are 
featured by a quantum oscillatory behavior. Furthermore, The prominent change in the Be(0001) 
surface electronic structure by the atomic hydrogen adsorption has also been shown. 

PACS numbers: 73.61.-r, 73.20.At, 73.21.Ac, 
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I. INTRODUCTION 

When the thickness of thin films approaches the 
nanoscale, the oscillatory quantum size effects (QSE) 
associated with electronic confinement and interference 
will occur— due to the splitting of the energy-level 
spectrum into subbands normal to the plane of the 
films. It has been shown that a change in film thick- 
ness by just one atomic layer can result in property 
variations on the order of 1/N, where N is the thick- 
ness of the film in terms of monolayers (ML). The os- 
cillatory QSE have long been clearly observed in ul- 
trathin metal overlayers on metal substrates^. Re- 
cent systematic experimental and theoretical investiga- 
tion of the QSE has mainly been focused to Pb films 
deposited on Si(111^^^03^^^^ or 
Cu(lll)^ substrates. 

In this paper we present a detailed first-principles 
study of the electronic structure and adsorption energet- 
ics of Be free standing films. The QSE in other simple- 
metal films such as Al 2 ^ 2 ^, Mg 24 ) L ^ Rh 24 ) and 
Pb 20 i 28 have been reported in previous references. Al- 
though Be is also a simple metal, its electronic struc- 
ture is profound compared to its close neighbors in the 
Periodic Table. In particular, the electronic proper- 
ties of bulk Be are very distinct from those in Be sur- 
face. For instance, the electron density of states (DOS) 
in bulk Be has a very deep minimum at Fermi level 
(Ep), making this material nearly semiconductin g 29 i 30 
while in the Be(0001) and Be(1010) surfaces the DOS 
at Ep was found to be larger by a factor of 4-5 that 
makes these surfaces two-dimensional free-electron-like 
metals 3 i. Compared to Mg and Li which have been pre- 
viously studied 2 ^ 2 - 7 in QSE, the bonding energy of Be is 
surprisingly strong (cohesive energy Be: 3.32 eV/atom, 
Mg: 1.51 eV/atom, Li: 1.63 eV/atom 3 ^) . Also the Debye 
temperature is very high (Be: 1440 K, Mg: 400 K, Li: 344 
K-2 3 .). Thus quantum effects should be more pronounced 
in determing low-temperature structural properties 3 ^. 

This paper is organized as follows: In Sec. II, the ab 
initio based method and computational details is out- 
lined. In Sec. Ill, the surface properties of the Be(0001) 



films, including the electronic structure, surface energy, 
work function, and interlayer relaxation, as a function of 
the thickness of the films, are presented and discussed. 
Also the properties of adsorption of atomic hydrogen 
monolayer onto Be(0001) surface is discussed in detail 
by presenting the sensitivity of the adsorption energy 
and local DOS to the thickness of the Be(0001) films. 
Finally, Sec. IV contains a summary of the work and our 
conclusion. 



II. COMPUTATIONAL METHOD 



The calculations were carried out using the Vi- 
enna ab initio simulation package 3 ^ based on density- 
functional theory with ultrasoft pseudopotentials 3 ^ and 
plane waves. In the present film calculations, free- 
standing Be(0001) films in periodic slab geometries were 
employed. The periodic slabs are separated by a vac- 
uum region equal to 20 A. In all the calculations be- 
low, a surface (lxl) was employed for the supercell 
slab. The Brillouin-zone integration was performed us- 
ing Monkhorse-Pack scheme 37 with allxllxlfc- 
point grid, and the plane-wave energy cutoff was set 300 
eV. Furthermore, the generalized gradient approximation 
(GGA) with PW-91 exchange-correlation potential has 
been employed with all atomic configurations fully re- 
laxed. First the total energy of the bulk hep Be was 
calculated to obtain the bulk lattice constants. The cal- 
culated a- and c-lattice parameters are 2.272 Aand 3.544 
A, comparable with experimental 3 ^, values of 2.285 Aand 
3.585 A, respectively. The use of larger fc-point meshes 
did not alter these values significantly. A Fermi broad- 
ening of 0.1 eV was chosen to smear the occupation of 
the bands around Ep by a finite- T Fermi function and 
extrapolating to T = K. 
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FIG. 1: (Color online) GGA energy bands and density of 
electron states (right panel) of hep bulk Be. The dashed line 
denotes Fermi level. 



III. RESULTS AND DISCUSSION 
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A. Band structure 

We first studied the properties of electronic structures 
of Be(0001) films. As a first step, we present in Fig. 1 
the band structure and the DOS of bulk hep Be. One 
can see that the density of electronic states for bulk Be 
resembles somewhat that of a semiconductor since it has 
a minimum at the Fermi energy. This makes Be differ- 
ent from its close neighbors, such as Mg, whose band 
property is nearly free-electron like. The bulk Be bands 
display a direct gap in a large range of the Brillouin zone, 
unlike Mg. Mg has a filled state at T with energy ~ 1.3 
eV— , while the corresponding Be state is above the Fermi 
energy and its band is nearly flat. This band is the source 
of both the low density of states near the Fermi energy 
and the high peak above the Fermi energy. Although the 
electronic configuration of elemental Be is ls 2 2s 2 , one 
can see from Fig. 1 that the p-orbital component in bulk 
Be plays a main role around Ep. The bonding properties 
of hep Be is anisotropic, which can be seen by the fact 
that the c/a ratio (1.56) is one of the most contracted 
for hep metals (for Mg, c/a ~ 1.62). Thus out-of-plane 
neighbors have shorter bonds than in-plane neighbors. 
Another evidence (not shown here) of this anisotropy is 
that the contribution of p x and p y orbitals to the DOS is 
different from that of p z component^. 

The above-mentioned "semiconducting" metallic pic- 
ture for bulk Be are totally changed in the case of 
Be(0001) surface. To illustrate this change, the elec- 
tronic structure properties of Be(0001) film of 10 layers 
are shown in Fig. 2, wherein the left panels give the re- 
sults for the unrelaxed structure, while the right panels 
show the relaxed results. For a film with the thickness 
of 10 layers, it shows in Fig. 2 no prominent difference 
in electronic structures between relaxed and unrelaxed 



FIG. 2: (Color online) Surface band structures and density 
of states of a Be(0001) slab with 9 Be monolayers. The left 
and right panels are for the unrelaxed and relaxed geometries, 
respectively. 



slabs. Compared to the bulk results shown in Fig. 1, it 
reveals in Fig. 2 that (i) the Fermi level in the Be(0001) 
film shifts down towards the broad peak, resulting in a 
prominent enhancement of the DOS at Ep with respect 
to the bulk case. This phenomenon of high density of 
surface states at Ep was previously reported^, and was 
also noticed in other metal films such as W(110) and 
Mo(110) 40 ; (ii) the band structure of the Be(0001) slab 
is characterized by a series of subbands which can be 
well fit by E n + h 2 (k1 + /c 2 )/2to* with the effective mass 
to* . This clearly demonstrates the quasi-2D free-electron 
character in the Be(0001) thin film, which is contrary to 
the bulk Be. 



B. Energetics 

Figure 3(a) shows the total energy per monolayer 
E(N)/N as a function of the number N of layers in the 
slab. The atoms in the slabs have been fully relaxed dur- 
ing the calculations. One can see from Fig. 3(a) that with 
increasing the thickness, E(N)/N gradually approaches 
a constant value which in the limit is equal to the energy 
per atom in the bulk Be. 

A quantity more suitably tailored to the QSE is the 
energy difference AE(N) = E(N) - E(N - 1). This 
quantity is shown in Fig. 3(b), which reveals damped 
oscillations as a function of N. These oscillations arise 
from the occupation of electronic levels close to Fermi 
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FIG. 3: (a) Monolayer energies E(N)/N for fully relaxed 
Be(0001) slabs as a function of the number of Be monolayers 
N; (b) corresponding energy differences . 



TABLE I: Interlayer relaxations given in percent, Adj,i+i, of 
Be (0001) as a function of the number N of layers in the slab. 



AT 


Arl 


A rl 


/\rl 


Arl 


AW 


A rl 


o 
2 


+ 10.9 1 2 












Q 



+6.540 


i a CAd 
+6.o4o 










4 


+5.4 


+2.2»8 


+0.4 








r 




1 Q 1 07 

-\-o.12i 


+U. I 00 


+U. 1 DO 


1 9 1 07 









+ Z.OOO 


fl 1Q1 

-U.lol 


+U. 1 02 


A 1 Q1 

-U.lol 


+ 2.000 




■7 


+z.4iy 


A Odd 

-0.206 


+U.U62 


+U.U62 


n nee 
-0.266 


+2.419 


8 


+1.718 


-0.767 


+0.176 


-0.29 


+0.175 


-0.767 


9 


+1.5 


-0.827 


-0.081 


+0.052 


+0.052 


-0.081 


10 


+1.032 


-0.829 


-0.001 


-0.025 


+0.173 


-0.025 


11 


+1.158 


-0.920 


-0.071 


+0.028 


+0.092 


+0.092 


12 


+0.935 


-1.008 


-0.003 


-0.214 


+0.277 


-0.073 


13 


+1.228 


-0.902 


-0.09 


-0.011 


+0.091 


+0.031 


14 


+1.038 


-0.951 


-0.031 


-0.21 


+0.216 


-0.073 


15 


+1.228 


-1.024 


+0.022 


-0.225 


+0.1789 


-0.018 



5.7 



u 5.4 



o 

tj 5.1 
c 



4.8 



4.5 





-a— Before relaxation 
^— After relaxation 



2 4 



10 12 14 16 18 20 22 24 
N 



FIG. 4: Work function of Be(0001) thin film as a function of 
the number of Be monolayers in the slab. 



surface?-, which are p-like as shown in Fig. 2. 

Figure 4 shows the work function of the Be(0001) thin 
films as a function of the number of Be layers in the slab 
for the unrelaxed and the relaxed geometries. One can 
see that the large influence of the relaxation on the work 
function occurs in the range of N ~ 3-7. For more larger 
values of N, the overall effect of geometry relaxation on 
the work function trend is small. Note that unlike the 
other metals such as Pb(lll), in which the work function 
oscillates periodically even in the ultra-thin film limit, 
the work function of Bc(0001) first curves up at small 
values of N (~ 3), then lowers down towards the bulk 
value in an oscillatory way. The peculiar behavior of work 
function in Fig. 4 implies the more complicated thin-film 
properties of Be compared to the other simple s-p metals. 
The similar behavior also occurs in the adsorption energy 
of atomic hydrogen onto Be(0001) surface (see below). 



C. Interlayer relaxation 



The interlayer relaxation, AcZ^+i, are given in percent 
with respect to the unrelaxed interlayer spacings, do, i.e., 
Adj^+i = 100(di.i + i — cIq)/cIq. di t i+i is the interlayer dis- 
tance between two adjacent layers parallel to the surface 
calculated by total energy minimization, do = c /2 is the 
bulk interlayer distance along (0001) direction. As men- 
tioned above, all layers in the slab were allowed to relax. 
Obviously, the signs + and — of Ad, v j +1 indicate expan- 
sion and contraction of the interlayer spacings, respec- 
tively. The relaxation of Be(0001) surface as a function 
of the number of layers is summarized in Table 1. Further- 
more, the interlayer relaxations are also plotted in Fig. 5 
as a function of N for clear illustration. One can see: (i) 
The three outmost layers relax significantly from the bulk 
value, in agreement with experimental observation^. In 
the whole range of layers that we considered, the topmost 
interlayer relaxation is always outward (Adi2 > 0), with 
Ag?i2 starting from the largest value of +11% for a slab 
with only two monolayers, and approaches a final value 
of ^1% with increasing the thickness of Be(0001) films. 
Note that the first interlayer separation on most metal 
surfaces is contracted, Be(0001) is one of the few excep- 
tions. Considering that the calculation is given at zero 
temperature, our calculated result of Ag?i2 ~ +1% can 
be considered comparable with the experimental observa- 
tion of Arfia ~ +3% at T=110 K. (ii) When the number 
of monolayers in the slab reaches N ~ 5, the interlayer 
spacings begin to oscillate with a damped magnitude. 
Again, the oscillations reflect QSE in Be(0001) thin films. 
After 25 monolayers, which is the maximal layers consid- 
ered here, the oscillations are still visible which indicates 
that the semi-infinite surface limit is still not reached. 
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FIG. 5: (Color online) Relaxations of the Be(0001) surfaces 
as a function of the number of Be monolayers in the slab 
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FIG. 6: (Color online) Adsorption sites for H adatoms on the 
hep Be(0001) surface. The indicated adsorption sites are on- 
top, hep, bridge, and fee. The H adatoms and Be atoms in 
the topmost surface layer are indicated by large open circles 
and small black circles, respectively, while the second-layer 
Be atoms are depicted as small red circles. 



D. Adsorption of atomic hydrogen: QSE of binding 
energy 

To further illustrate the physical properties influenced 
by finite size of the thin films, in this section we focus 
our attention to the adsorption of atomic hydrogen on 
Be(0001) thin films. To the best of our knowledge, the 
reflection of QSE by the adsorption features has not been 
previously studied. 

Before we study the hydrogen adsorption properties 
as a function of the thickness of Bc(0001) thin films, we 
need to determine the energetically favourable adsorp- 
tion site. Since the preference of adsorption site is not 
sensitive to the thickness of the substrate, thus to look 
for this preference of the adsorption site, it is sufficient 
to give a study on the slabs with fixed thickness of the 
Be(0001) substrate, which at present is chosen to be 9 
monolayers. We choose four most probable adsorption 
sites, namely, on-top, bridge, fee, and hep sites, which 
are indicated in Fig. 6. The binding energy is calculated 
using the following equation: Binding energy [atomic 
H]= -(£[H/Be(0001)] - B[Be(Q001)] - 2E[K])/2 where 
_E[H/Be(0001)] is the total energy of a slab which con- 
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FIG. 7: Calculated binding energy of H adatom as a function 
the number of Be(0001) monolayers in the slab. 



sists of 9 layers of Be atoms and one H atom on each side 
keeping inversion symmetry, 2?[Be(Q0Ql)] is total energy 
of the slab without H atoms, and E[H] is total energy of a 
free H atom which is put in a 10 Ax 10 Ax 10 Asupercell. 
As a result, the calculated binding energy of atomic H for 
different adsorption configurations is 1.46 eV (on-top), 
2.24 eV (fee), 2.11 eV (hep), and 2.33 eV (bridge), show- 
ing a clear preference for bridge site adsorption. We note 
that this bridge-site preference for H adsorption was also 
reported in Ref. [43] and in other metal surfaces such as 
W(100)ii. 

After finding the preferred atomic H adsorption site 
(bridge site), we give a series of calculations for the bind- 
ing energy of the H adsorbate as a function of the thick- 
ness of Be(0001) thin films. The results are summarized 
in Fig. 7. One can see that the binding energy curves up 
at small film thickness, followed by damped oscillations 
when increasing the Be monolayers in the slab. Thus the 
binding energy of atomic H depends on the thickness of 
the quantum films in an oscillatory way. In experiment 
this dependence can be observed by investigating the de- 
pendence of H coverage on the monolayers of Be(0001) 
thin films. Note that the dependence of H adatom bind- 
ing energy on the thickness of the Be(0001) films looks in 
analogy with the work function behavior shown in Fig. 
4. The cause of this similarity remains unclear at the 
present stage. 

The oscillatory thickness dependence of the binding 
energy should be related with the bonding properties be- 
tween H adatom and surface Be atoms. To obtain a 
deeper understanding of the physics behind this obser- 
vation, we turn to a full study of the surface electronic 
structure of H-covered Be(0001). For this we calculate 
the local density of states (LDOS) through projections of 
the total wave function onto atoms of interest within the 
Wigner-Seitz spheres around them. The Wigner-Seitz 
radii for the H and Be atoms are taken to be 0.5 Aand 
1.0 A, respectively. Figure 8 plots a series of LDOS of H 
adatom layer and the topmost Be(0001) layer for differ- 
ent thickness of Be thin films. A strong chemical bonding 
of characteristic molecular H-Be bond can be seen. In 
the case of ultra-thin slab with only one Be monolayer 
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FIG. 8: (Color online) Calculated local density of states of 
H adatom and topmost Be atom. The number of Be(0001) 
monolayers in the slabs is shown in the figure. 



contained, one can see that the Fermi level falls in the 
broad-peak region of Be surface. Since the Fermi level for 
the clean Be(0001) surface also falls in this broad-peak 
location (see Fig. 2), thus it suggests that adsorption of 
atomic H onto Be monolayer does not change the surface 
DOS much. As a consequence, the H binding energy is 
lowest in the case that only one Be monolayer is included 
in the slab (see Fig. 7). When the Be(0001) monolayers 
in the slab is increased, one can see from Fig. 8 that the 
Fermi energy continuously shifts upwards, approaching 
a narrow peak, showing a more strong hybridization be- 
tween H Is and Be 2p orbitals. Consequently, the binding 
energy of H adatom first curves up by increasing the Be 
layers in the slab, as shown in Fig. 7. When the thickness 
of slab increases up to 5 Be monolayers, then the Fermi 
level remains almost unchanged with respect to the nar- 
row peak, the resultant hybridization between H Is and 
Be 2p orbitals obtains its largest amplitude. Therefore, 
when further increasing the layers of Be in slab, instead 
of a monotonic increase or decrease, the H adatom bind- 



ing energy is now featured by an oscillatory behavior due 
to metallic nature of Bc(0001) surface. Furthermore, by 
a comparison of the clean-surface DOS (Fig. 2) with the 
H-adsorbed surface DOS (Fig. 8), remarkably, one can 
see that there is a prominent change in the Be(0001) sur- 
face DOS by the atomic hydrogen adsorption. Since it 
has been speculated that the metallic enhancement of Be 
surface is responsible for many extraordinary phenomena 
such as enhanced electron-phonon coupling 4 ^, the 400 
increase in the superconducting temperature for amor- 
phous Be compared with crystalline Be, clearly it is of 
high interest to study the effect of H adsorption on these 
peculiar surface-related phenomena. 

IV. CONCLUSION 

In summary, the Be(0001) thin films have been studied 
by density-functional theory pseudopotential plane-wave 
calculations. The dependence of electronic structure, en- 
ergetics, and interlayer relaxation upon the thickness of 
the film has been fully investigated, showing the metallic 
QSE of the film, although the bulk Be is nearly semi- 
conducting. We have shown that this QSE is associated 
with an interplay of surface enhancement at the Fermi 
level and the bulk band splitting into subbands due to 
the quasi-2D electronic confinement. We have also stud- 
ied the atomic hydrogen adsorption onto the Be(0001) 
films. The results gave a clear preference for bridge site 
adsorption. We have shown that (i) the H adsorption en- 
ergy oscillates with increasing the Be monolayers in the 
slab; (ii) In the case of ultra-thin films, the H adsorbate 
changes the surface density of states prominently by con- 
tinuously shifting the Fermi level upwards, approaching 
the narrow peak in the otherwise clean-surface DOS. The 
shift will saturate after the thickness of Be(0001) film in- 
creases up to 10 layers. We expect this prominent change 
in the density of states will have a significant effect on the 
other relevant physical properties in Be(0001) surface. 
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